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We analyze in detail the analytical solutions of the Dirac equation with scalar S and vector V 
Coulomb radial potentials near the limit of spin and pseudospin symmetries, i.e., when those po- 
tentials have the same magnitude and either the same sign or opposite signs, respectively. By 
performing an expansion of the relevant coefficients we also assess the perturbative nature of both 
symmetries and their relations the (pseudo)spin-orbit coupling. The former analysis is made for 
both positive and negative energy solutions and we reproduce the relations between spin and pseu- 
dospin symmetries found before for nuclear mean-field potentials. We discuss the node structure 
of the radial functions and the quantum numbers of the solutions when there is spin or pseudospin 
symmetry, which we find to be similar to the well-known solutions of hydrogenic atoms. 
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I. INTRODUCTION 

Spin and pseudospin symmetries are SU (2) symmetries of a Dirac Hamiltonian with vector and scalar potentials. 
They are realized when the difference, A = V — S, or the sum, £ = V + S, are constants. These constants are zero for 
bound systems whose potentials go to zero at infinity. Generally, in relativistic quantum physical systems with this 
kind of potentials neither of these conditions is met exactly but, in some cases, one of them can be approximately true. 
These symmetries may explain degeneracies in some heavy meson spectra (spin symmetry) or in single-particle energy 
levels in nuclei (pseudospin symmetry), when these physical systems are described by relativistic mean- field theories 
with scalar and vector potentials However, for bound systems whose potentials go to zero at infinity, pseudospin 
symmetry cannot be realized, because £ = S + V (which must go to zero in the limit of exact pseudospin symmetry) 
is also the binding potential. For anti-fermions the opposite occurs, i.e., one cannot have exact spin symmetry, i.e., 
V = S 0, H. This behavior has been related by several authors to the perturbative nature of spin and pseudospin 
symmetry [J-Q- 

The problem of a fermion moving in the background of a general mixing of vector and scalar Coulomb fields has 
been originally solved for investigating the implications of a tiny contribution of the scalar Coulomb potential to the 
atomic spectroscopy @ . In this paper we explore the fact that this problem has analytical solutions for arbitrary (to 
a point) scalar and vector Coulomb radial potentials to examine in a detailed and explicit way, both for fermions and 
anti-fermions, the onset of spin and pseudospin symmetries and assess their perturbative nature. 

This paper is organized as follows. In Sec. [TT| we review the general solutions of the Dirac equation with spherical 
scalar and vector Coulomb potentials, discuss their quantum numbers and their relation with (pseudo)spin quantum 
numbers, and present the formulas for the energy eigenvalues. We also discuss in detail the node structure of the 
radial wave functions and compare them to what is known from other kind of radial mean-field potentials. In the 
next section we present the expansions of the energy eigenvalues for positive and negative energy solutions in terms 
of relevant potential parameters for spin and pseudospin symmetries and discuss the perturbative nature of those 
symmetries. Next we discuss the quantum numbers and extra degeneracies of the Dirac equation solutions when there 
is spin and pseudospin symmetry. We find that these are similar to the ones found in non-relativistic hydrogenic 
atoms. Finally, in Section HVl we draw the conclusions. 

II. SOLUTIONS OF THE DIRAC HAMILTONIAN WITH SCALAR AND VECTOR COULOMB 

POTENTIALS 

The Dirac Hamiltonian with scalar S and vector V potentials reads 



H = a ■ pc + f3(mc 2 + S) + V , 



(1) 
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where a and /3 are the Dirac matrices in the usual representation 

er\ 
.0 ' 



I 
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where <x are the Pauli matrices and / is the 2x2 unit matrix. If S and V are of Coulomb type one has 

ay 



V = 



he 



S = ^he 
r 

The Hamiltonian ([!]) can be written in terms of the sum and difference potentials £ = V + S and A = V — S as 



H = apc + /3mc 2 + -(I + /3)E + -(I - (3)A 



where 



£ = — — He 



r 



he 



«s = av + as 
aA = ay — as 



The solution of the time-independent Dirac equation 



with potentials © and © is of the form 



Hi, = Eip , 
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where 



-(1+1) j = l+\ 
I 3=i-\ 



(10) 



I is the orbital angular momentum of the upper component and k = —re. The angular functions <p Kmj {8,ip) are the 
spinor spherical harmonics and g nK (r) and f K (r) are the radial wave functions for the upper and lower components of 
the Dirac spinor respectively. The orbital and total angular momenta can be obtained from re by I = |re| + l/2(re/|re| — l) 
and j = \k\ - 1/2. 

The analytical bound-state solutions of eq. (JSJ with the Hamiltonian (J5]) with the potentials (O and (JTJ) can be 
taken from Greiner and Rafelski book [7j|, replacing the scalar and vector potentials by their sum and difference. One 
has for the radial wave functions, using the notation of Leviatan [8] 



9n(r) 
Mr) 



-A^/mc 2 + E [(k + m)Fi+ n r F 2 ] p 1 e~ p/2 



Ay/ me 2 - E [(k + m) Fi ~ n r F 2 ] p 1 e~ p/2 
where F\ = iFi(—n ri 2"/ + 1, p) and F 2 = \Fi(—n r + 1, 27 + 1, p) are confluent hypergeometric functions and 
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The integers n r — 0,1,2,... are the quantum numbers defined by the quantization condition (|16p required in order 
that the radial functions g K (r) and f K (r) be normalizable, and the respective normalization condition 



L 



oo 

[gl(r) + fl(r)]dr = l. (19) 
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determines the constant A. 

One usually defines the principal quantum number as n = n r + \k\ n = 1, 2, ... so that one can write (|16p as 

»fc = -n+|/s|-7=-£. (20) 

Since £ = n r + 7 is positive definite, this last relation, together with the condition that A, given by (|14l) . must be real 
to have bound solutions, yields the following constraints to the energies and coefficients a a and as: 

\E\ < mc 2 (21) 
a^{E + mc 2 ) + a A (E - mc 2 ) < . (22) 

This last equation can be rewritten as 

mc 2 — E 

as < a A r — . (23) 

mc z + E 

Another condition comes from the requirement that 7 be real. One has 

a E a A < k 2 . (24) 



This condition and (|23[) imply that there is no bound-state solution if «s > and a A < 0. If as and a A have the 
same sign both conditions (|23l) and (|24l) constrain the values of the strengths of the potentials. When as < and 
aA > those conditions do not restrict the strength of these potentials. On the other, since one would also like to 
have ay < (standard attractive Coulomb potential) for positive energy states, then one must have |as| > |«a| if 
the potentials have diferent signs. Thus, in the examples presented below, we shall use values for the strengths of the 
Coulomb potentials such that as < 0, a a > and — as > ctA- 

From (fH)|) . (fTT|) and (|T^|) and one can derive the following useful relations: 

vl = Vi+ "Atts (25) 
— n r (n r + 27) + k 2 

From eqs. (|T4|. (fT7|) and (j20|) one can calculate the eigenenergies. One gets two types of solutions, denoted by E^: 



Et R = mc 2 al-^^^Ve+p az 

(a A +as) 2 +4£ 2 v ; 



The dependence on n r and k comes through £ = n r + 7 = n r + \/ k 2 — oa»s = n — \k\ + V^ 2 «a«s- 

The signs were chosen such that i? + is a positive energy state and E~ is a negative energy state when either cka or 
as are very small (see the following section). Since charge conjugation ("c.c") produces the changes V —> —V and 
S ->• S, E ->• -E 1 , one has 

c c c c c c 

aA '> —as as '> — a a £ "> £ (same quantum numbers k and n r ) 

+ cc " W 

This means that the whole spectrum maps into itself with respect to the charge conjugation operation and therefore 
contains both particle and anti-particle states. Also we note that there is a double degeneracy in the energy levels 
(besides the spin degeneracy) , because the eigenenergies (|26[) only depend on the magnitude of k (except when n r = 
— see the discussion below). As remarked by Leviatan [8] this degeneracy has to do with the fact that the Hamiltonian 
([1]), when the scalar and vector potentials are of Coulomb type, commutes with the operator 

S • r 

J = -U<~f 5 (H ~ I3mc 2 ) + (a v mc 2 + a s H) (28) 

r 

= -iA7 5 (i?-^mc 2 ) + ^-^[(as + aA)mc 2 + (as-a A )H]/2 . (29) 

r 



4 



where 

a ) (30) 

and K is related to the L S operator (see section lPL 7 )) such that the spinor ^ is its eigenspinor with eigenvalue n. This 
is a generalization of the Johnson-Lippmann operator [8, 9] , which is itself a generalization of the Laplace- Runge-Lenz 
operator of the non-relativistic quantum Coulomb problem to the relativistic one. Since J anti-commutes with K , 
the energy must depend only on 

One may check that the non-relativistic limit is correct. If one sets ag = «a = &v an d have ay -C £ , £ ~ n, one 
gets 

2 

E + -mc 2 mc 2 ?% (31) 

2n 2 

Given that, from ([3]), we may set, for the quantum system of an electron in a atom of atomic number Z, ay = —Za, 
where a = e 2 /(hc) is the fine structure constant, equation pip is just the energy of the level n of an electron of a 
hydrogenic atom with atomic number Z. 



A. Node structure of the solutions 



We discuss now the node structure of the radial wave functions (fTTj) and (|12p. The solution with n r = deserves a 
special consideration. In this case the radial functions g K (r) and f K (r) will reduce to 



g K (r) = -B\Jmc? + E{k + r/ 2 ) p 1 e~ p/2 (32) 
f K (r) = Byjmc 2 - E (k + r] 2 ) p* e - p/2 (33) 

where B is a constant. These functions have only one node at the origin. Furthermore, we may note that, if as and 
Qi have the same sign, 7 < 1, and therefore the radial wave functions g K (r)/r f K (r)/r, although integrable, would 
be singular at the origin. Another important point is that, in this case, from (1251) . 772 = ±H- Because of the factor 
K + r/2 in eqs. (I32[) and (1331) a non-zero wave function cannot exist for any sign of n. It turns out, as already remarked 
by Leviatan [8|, that for E + solutions one must have k < while for E~ solutions one must have k > 0, which means 
that for n r = the degeneracy in the sign of k is broken. 
For n r > 1, g K (r) and f K (r) can be written as 



g K (r) = -AWmc* +E[(k + m ) I%{p) + (2 7 + n^L^ip)] e"^ 2 (34) 
f K (r) = A'Vmc^-E [(« + r, 2 )L 2 r7M - (27 + «r)^_i(p)] P 1 e~ p/2 (35) 

where L 2 J(p) are the generalized Laguerre polynomials of degree n r = n— \k\ and A' is a constant. These polynomials 
have n r distinct positive zeros. These formulas suggest, then, that g K (r) and f K (r) may have the same number of 
nodes, namely n r + 1 (n r = n — \k\ nodes for r > plus one at r = 0). In the following we prove that, if a^aA < 
then indeed the number of nodes of g K (r), n g , and the number of nodes of f K (r), n/, are both equal to n r + 1 for any 
sign of k. 

The zeros of generalized Laguerre polynomials of different degrees are interweaved, i.e., the ith zero of L 2 ^ _ 1 (p) 
is localized between the ith and (i + l)th zero of L 2 J(p) and their value at p = is positive and increases with the 
polynomial degree. As a consequence, one finds that the number of zeros riz, of the linear combination (A / 0, 5 ^ 0) 

L(p)^AL 2 Z(p) + BL 2 n l_ 1 (p) (36) 

depend on the ratio B/A in the following way 

B ^ J#(0) 
V - 1 , -r < 



5 > ^(o^ • (37) 

A ~ inl-l(O) 
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From the values at the origin of associated Laguerre polynomials (see, for instance, [T(|) one has 

L%(0) _ n r + 2 1 



(38) 



Since there is an extra node at r = p = and that the equality in f|3T[) corresponds to a node at the origin, from 
eqs. and (|3"5|) one may write 

27 + n r n r + 27 
V , ; < 

27 + n r n r + 27 

lr + 1 , > 

K + i] 2 n r 
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sgn(« + 772)n r 


n r -+■ 
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sgn(« + r] 2 )n r 


n r 




|k-1 






m(K + 77 2 )n r 


n r -+■ 


•1 , 


|/« 4 


- »72 1 
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m(K + rj 2 )n r 


I the 


sign 


of X. 









This can be further simplified to get 

[ n r , \k + rj 2 \ <— sgn(n + rj 2 )n r 

(39) 

n f = i , (40) 

[ n r + 1 , |k + 772 1 > sgn(K + r] 2 )n r 

tion giving the sign of x. 
On the other hand, from ([T5|) and (|23|) and when as a a < 0, the following inequalities hold 

\r] 2 \ > n r 

7 > |«| 1 772 1 > "r + |«| 

meaning that 

|/t + 772|>JV (41) 

which is true when k and 772 have the either the same or different signs, which is to say that it is true for any sign of 
k and 772. From (f3"9"]l and (|4"0"1) it follows that, when as«A < 0, one has 

n g = 7i,. + 1 (42) 
nf = 77,. + 1 (43) 

for any k. 

Note that this result applies to both E + and E~ solutions. Also, since as and «a have different signs, 7 > 1 and 
thus g K (r)/r and f K {r)/r have the same nodes as g K (r) and /«(r). 

In Fig. [T]this is shown explicitly for the E + levels 2si/2, %Pi/ 2 (k = =pl, n r = 1) (left panel) and for the E~ levels 
4^3/2, 4d 3 / 2 (k = T2, n r — 2) (right panel), for oa = 0.5 and as = —0.8. The radial coordinate is given in units 
of the Compton wavelength Lq = h/(mc). Here one uses, as usual, the spectroscopic notation n£j referring to the 
upper component of the Dirac spinor. 

It is interesting to remark that this node structure does not follow the rule for the Dirac equation with vector and 
scalar radial potentials derived by Leviatan and Ginocchio in [Tl| . 

{7l„ , k < 

(44) 
n g + 1 , K > 

by which only for negative k the upper radial and lower radial functions of the Dirac spinor have the same number 
of nodes. The structure (pHf also holds for anti-particles in Woods-Saxon scalar and vector mean-field potentials, as 
shown in 3j- This may be because these potentials were finite at r = 0, as are the nuclear mean- field scalar and 
vector potentials, and thus the behavior of the respective radial functions at the origin is different form the present 
case of Coulomb potentials. 
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Figure 1: Radial wave functions <?n(7')/r and f K {r)/r for the E + states 2si/2, %Pi/2 ( K = Tb n r = 1, n g j = 2) and the E 
states 4p 3 / 2 , 4d 3 / 2 (« = T2, TV = 2, n 9 j = 3), when qa = 0.5 and as = —0.8. he = h/(mc) is the Compton wavelength. 



III. SPIN AND PSEUDOSPIN SYMMETRIES IN THE DIRAC HAMILTONIAN WITH SCALAR AND 

VECTOR COULOMB POTENTIALS 



To study the conditions under which the Dirac Hamiltonian ([5} with the Coulomb potentials (JBJ and ([7]) has bound 
solutions near the conditions of spin and pseudospin symmetry, we will perform an expansion of eq. (|26[) for small 
a a (spin) and as (pseudospin). 

If we expand (|26[) (using the principal quantum number n) until next-to-next-to-leading order in oa one gets 



mc < 1 



24 



4a| (k 2 - 2n|«|) 
"4^ K 2 (q| + 4n 2 ) 2 



aA 



t S 1 rv ^u^ 

a| {4n 3 | K | [a| + 4 (n 2 - 4 K 2 )] + k 2 [a 2 {k 2 - An 2 



20n 2 n 2 + A8n 4 ] } 



2n A n 2 (a 2 + 4n 2 )- 



E„ „ = mc 



-1 



a\ a s (k 2 - 2n\n\) 



2n 2 



An 2 n 4 



A) 



while for as one has 



E+, 



a| aA (k 2 — 2n\n\) 



2n 2 



A K 2 n 4 



a 



e: 



mc < — 1 



+ - 



2ai 



"A 



4n 2 



k 2 (a 2 ^ + 4n 2 ) 



1 + 0(4) 

2n|«|) 
2- 



2k% 2 (a 2 , + 4n 2 ) 3 



a 2 A +0 (a A ) 



{4n 3 | K | [a 2 ^ + 4 (n 2 - 4 K 2 )] + k 2 [a 2 , ( K 2 - 4n 2 ) + 2fj K 2 n 2 + 48n 4 ] } 



•O (a 




(45) 
(46) 

(47) 



(48) 



From these equations is clear that E + ^ are the energies of bound positive (negative) energy states for small a a and 
as- The symmetry of the expressions above when one switches from aA to as and from positive to negative energy 
states is just the illustration of the charge conjugation transformations as given by eqs. (|27p . From the equations 
above one sees also that when the zeroth-order term is just ±mc 2 the respective symmetry cannot be realized, since 
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there is no bound states. For positive energy states one sees that in eq. (14"7T) (pseudo-spin symmetry) and for negative 
energy states in eq. (|46[) (spin symmetry). Again, this agrees with what is known from realistic calculations in nuclei, 
since E and —A potentials act as biding potentials for fermions and antifermions respectively @, Q. 

It is interesting to examine how the spin and pseudospin symmetries get broken and the corresponding energy 
splittings of the (pseudo)spin partners. 

The next to leading order terms of equations (|45M8p reveal that the breaking of spin symmetry is perturbative 
and the breaking of pseudosymmetry is non-perturbative for E + states. The reverse is true for E~ states. Indeed, 
this can be seen the fact that the next to leading order terms in eqs. (|4"5j) and (|47l) are first-order in a a and as, 
respectively. This is what we expect from first-order perturbation theory and also from the fact that the breaking of 
spin and pseudo-spin symmetries from their exact realizations should come from the spin-orbit and pseudospin-orbit 
coupling terms, which can be obtained from the second-order radial equations for g K (r) and f K (r), respectively (see, 
for instance, ref. Q) 



(he) 



(he)' 



1 d 
r 2 dr V dr 



1 d 

r 2 dr 



dr 



k(k + 1) 



k(k — 1) 



A' 



E + mc 2 — A(r) \dr 



mc* 



E(r) \dr 



1 + K 



9n(r) 



[E - mc 2 - T,(r)][E + mc 2 - A(r)} 



1 — K 



-IE- 



f«(r) 



E(r)P- 



A(r)] 



9n(r) 



Mr) 



(49) 



(50) 



where the prime means derivative with respect to r. These are proportional to the derivatives of A (spin) and E 
(pseudospin) and thus, for small enough values of a a and as, respectively, their expectation values contributing to 
the energy splitting of spin and pseudospin partners should be linear in those parameters. In particular, one may 
relate the contributions from spin-orbit and pseudo-orbit couplings to the total energy to the following terms coming 
coming from eqs. (|4T)]) and ([5TJ)) as was done in ref. 



o 







These terms are consistent with the non-relativistic formula for spin-orbit coupling, the first one applied to positive 
energy states and the second one to negative energy states. The origin of the factors 1 + n for the upper component 
and to 1 — k for the lower component can be traced to the action of the operator L ■ S upon the spinor since 
L ■ Sijj — —h 2 /2(/3K + I)tp with Kip = nib and S = h/2 E. However, contrary to what happens with Woods-Saxon 
type potentials for positive energy states the spin-orbit term ((STj) has a double pole at r = hca^j(E + mc 2 ) so it 
cannot be calculated separately from the so-called Darwin term, coming from the derivative term in eq. (|49[) and also 
multiplied to A', since the sum of the two terms must of course be finite. 

Therefore, for vector and scalar Coulomb potentials, one can establish a connection between the realization of exact 
pseudospin and spin symmetries for bound states and their perturbative breaking. This was also recently shown to 
be the case for vector and scalar Woods-Saxon nuclear mean-field potentials for neutrons and anti- neutrons @. 



A. Quantum numbers for exact spin and pseudospin symmetries 

Finally, we discuss the quantum numbers of states when there is either exact spin or pseudospin symmetries, i.e., 
when oa = or as — respectively. The energies of those levels are given by 

< K =- 2 (l-^,) (spin) (53) 
En, K = mc 2 (-l + a 2 2 °l\ n ^ ) (pseudospin) (54) 
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One may notice immediately that these energies only depend on the principal quantum number n. This situation 
looks very similar to what happens with the energy levels of a non-relativistic hydrogenic atom. To see how it comes 
about we start out with spin symmetry case (A = or a a = 0). The corresponding Dirac equation reads 

[a ■ p c + f3mc 2 + ^(1 + = Ei> , (55) 

with S = a^hc/r. 

Following the procedure of [l2| we use the projectors ip± = P±ip = [(I ± /3)/2] ip such that 

* + = (°) *- = (")' <56) 

where <j> and \ are respectively the upper and lower two-component spinors. Using the properties and anti- 
commutation relations of the matrices j3 and a. we can apply the projectors P± to the Dirac equation (|55|) and 
decompose it into two coupled equations for tp + and ip-\ 

ca.-p4i_ + (mc 2 + S)-0+ = Ei/>+ (57) 
col ■ ptp+ — mc 2 ip_ = Etl'- ■ (58) 

From this last equation we get and replace it in (|57p to give 

p 2 c 2 i> + = {E + mc 2 ){E-mc 2 -Y>)4> + . (59) 

or yet 

|-^ + = (f'-S')^+, (60) 

where 

f - , L . _ 

\2mc 2 J ' \2mc 2 

Equation (|60p is just the Schroedinger equation for a hydrogenic atom of "atomic number" 

Z' = -(-^ + l)ax/a , (62) 
where again a is the fine structure constant. The energy of this "atom" is just (see eq. ([31]) ) 

Replacing the expressions (ICT]) into this last equation one gets eq. (|53|) . 

The fact that the energy only depends on n in spin symmetry conditions implies also that there is an extra level 
degeneracy. Because of the relation n = n r + \k\, levels with the same principal quantum number and \k\ < n are 
degenerate. This also has a correspondance to the well known degeneracy of levels of a hydrogenic atom in which all 
levels with £ < n — 1 and the same n are degenerate. These properties can be traced back to the fact that in spin 
symmetry conditions and radial potentials there is an SU(2) symmetry generated by the operator 



£'=fi + lV, S'=f-^ + l)S, £=E-mc 2 . 








U P LU P 



(64) 



where U p = <r ■ p/(\/p 2 ) is the helicity operator. This means that the orbital angular momentum £ of the upper 
component of the Dirac spinor is a good quantum number is this case 

mm. An ° ther 

way of stating this is saying 

that when there is spin symmetry there is no spin-orbit coupling. 

A similar reasoning can be made when there is pseudospin symmetry (E = 0, ). In this case one would get the 
following equation for the lower component : 

fU_ = (£'-A')V- , (65) 
2m 
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where 

2mc 2 / \ 2mc 2 

As one would expect for negative energy states and from charge conjugation, this amounts basically to change the 
energy sign and also replacing £ by —A [3j. The SU(2) generator would be in this case 

- 1 / U p L U p \ 

C= -r a-pLa-pP++LP- = , 67 

V \ L J 

and in this case the orbital angular momentum I of the lower component of the Dirac spinor is a good quantum 
number [H. \v\. 



IV. CONCLUSIONS 



We analyzed in detail the analytical solutions of Dirac equation with scalar S and vector V Coulomb radial potentials 
in terms of the sum £ and difference A potentials. Besides reviewing the solutions at this light, we established the node 
structure of their radial solutions and examined the solutions near the limit of spin and pseudospin symmetries, i.e., 
when when A = and £ = respectively. We assessed the perturbative nature of those symmetries, and confirmed 
previous results with another type of potentials, namely, that spin symmetry is perturbative and pseudospin symmetry 
in non-perturbative for positive energy solutions, while the reverse happens for negative energy solutions. We also 
found that in conditions of exact spin and pseudospin symmetries the solutions of the Dirac equation have similar 
features of the well-known non-relativistic solutions for hydrogenic atoms, namely that they depend only on the 
principal quantum number and exhibit a similar degeneracy. 
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